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Abstract 

We study a class of quadratic, infinite-dimensional dynamical 
systems, inspired by models for viscoelastic fluids. We prove that 
these equations define a semi-flow on the cone of positive, essentially 
bounded functions. As time tends to infinity, the solutions tend to 
an equilibrium manifold in the L 2 -norm. Convergence to a particular 
function on the equilibrium manifold is only proved under additional 
^ assumptions. We discuss several possible generalizations. 
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1 Introduction 



This paper is concerned with evolution equations of the form 
dy 

± = yna-y), y(x,0) = y (x), (1) 
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where y(-, t) is an unknown and fl(-) is a given real-valued function, both 
defined on a measure space (Q, fi) with finite mass (/-i(Q) < °o). The operator 
7 is an orthogonal projection on the Hilbert space L 2 = L 2 (Q /( u). We will 
use the standard notations (•, •) for the inner product in L 2 , and || • || p for 
the U = U(Q., y.) norms. We denote by L°° ,+ the cone of strictly positive 
functions in L°° = L°°(Q, fi): 

L°°' + = {u € L°° : ess inf u(x) > 0}. 

Equation ([I]) is subject to the following assumptions: 
Assumption 1 

(i) The operator !P : L 2 — > L 2 is an orthogonal projection, satisfying 
J > (L°°) c L°°. Its null space, ?sf(CP), is one-dimensional, spanned by an 
essentially positive function n e L°° /+ , which we take to be normalized, 
(n,n) = 1. 

(ii) The function a(x) € L°°. Without loss of generality we can assume 
that 

9(a) = fl. 

The system (H) is a toy model inspired by models of viscoelastic fluids. 
Specifically the Maxwell constitutive equation for the conformation tensor 

is in 

+ (« • V)a = (Vufo + o(Vu) + -(o - I). (2) 

at A 

Here a(x, t) is the symmetric positive-definite conformation tensor, u(x, t) 
is the velocity field, I is the unit tensor and A is the elastic relaxation time. 
For polymers suspended in an incompressible solvent in the creeping flow 
regime, the velocity field is derived from the conformation tensor through 
the Stokes equations. The mapping a h» Vm, is linear and can be explicitly 
written by means of a Green function, 

Vw(x) = I G n (x, y) ■ diva(y) dy, 
Jq 

where Q is the domain (which may be bounded on not) and Gq is the 
corresponding Stokes kernel (i.e., the Green function of the Stokes problem). 
It can easily be shown that this mapping is, in fact, an orthogonal projection 
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(see 13). Since the terms in equation (0 that can potentially lead to finite- 
time blowup are the stretching terms, it is of interest to omit the advection 
and the relaxation terms (which by themselves are not norm increasing), 
and consider systems of the form 

= [Ma - o)] T o + oy(a - a). (3) 

at 

Here a is an external force field (see [2 J for more details). The system ((D is 
a one-dimensional scalar toy model, that mimics the dynamics (3). 

Equation Q} can also be viewed as an infinite-dimensional generalization of 
a Lotka-Volterra system [31. In Section [2] we prove that defines a global 
(in time) semi-flow on the cone of positive functions L°° ,+ (Theorem 12.21 in 
Section[2]). We then proceed to analyze the long-time behavior of this system. 
It is clear that every function y satisfying 1 > (y) = IP(a) is an equilibrium 
solution of (dJ, and these are the only equilibria in L°° /+ . Our main theorem 
asserts that the equilibrium set 

is the global attractor for all initial data yo e L°° ,+ (Theorem B.ll in Section|3|. 
The convergence of y(-, t) to the manifold M, as t — > oo, is in the L 2 -norm. 
The theorem does not guarantee uniform convergence, nor does it guaran- 
tee that y(-, t) converges to a specific equilibrium in M. For this to happen, 
additional assumptions are made; various situations are considered in Sec- 
tion|U We conclude this paper with a discussion about open questions and 
various generalizations. 

2 Global existence 

We start by establishing the well-posedness of equation {T) under Assump- 
tion [TJ The first step is to show existence and uniqueness of solutions for 
short times: 

Theorem 2.1 (Local-in-time existence and uniqueness) Letyo e L°° be given. 
Then there exist times T\, T2 > 0, depending on yo only, such that (d) has a unique 
solution yeC 1 ((-T 1/ T 2 ),L°°). 

Proof: Note first that due to Assumption [T] the operator 7 has the explicit 
form 

y(z) = z - (z,n)n. 
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It is a bounded linear operator L°° — > L°° since 

ITOIU < (1 + K")INli)INU. 

We rewrite |(T]) as 

^| = ya-y:P(y)EEF(y). 

The short-time existence and uniqueness of solutions follows from Picard's 
theorem over Banach spaces, provided that F is a locally Lipschitz continu- 
ous mapping L°° — > L°°. This is indeed the case as 7 is a bounded operator, 
hence it is locally Lipschitz, and the product of locally Lipschitz functions 
is again locally Lipschitz. ■ 

We then show that solutions that are initially positive remain so at all times: 

Proposition 2.1 (Positivity) Let y e C 1 ((-Ti, T2); L°°) be a solution of CD, with 
initial condition yo e L°° ,+ . Then y(x, t) remains positive, i.e., ess inf x€ Q y(x, t) > 
0, for all t € (— T\, T2). In other words, the cone L°° ,+ is an invariant set for the 
dynamics. 

Proof: The positivity follows readily from the fact that the unique solution 
of ([I]) solves the integral equation 

y(; t) = y exp 3>(« - y(-, s)) dsj . (4) 



The next step is to show that the solution with initial data in L°° ,+ , as long 
as it exists, is bounded, uniformly in time, in L°°, by a constant that only 
depends on the initial data. The proof relies on the fact that the dynamics 
([l]) subject to Assumption[TJpreserve the natural order among functions. To 
simplify notations, we define Q := J - T to be the orthogonal complement 
of the projection "P, namely, Qy = (n, y) n. 

Lemma 2.1 Let y e L°° be a non-negative function, y > 0. Then, 

essinf Qy(-) > 0, 



with equality if and only ify-0. 
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Proof: The non-negativity of y and the positivity of n implies that 
essinf Qy(-) = (n, y) ■ essinf n(-) > 0. 

Since ess inf n{x) > equality occurs if and only if (n, y) = 0, i.e., if and only 
if y = 0. ■ 

Proposition 2.2 (Comparison principle) Let y,z e C 1 ([0, T);L°°' + ) be two so- 
lutions of (Q} with initial data yo,zo e L°° /+ . If yo > zo zn ^ 

y(v0>z(',0 (5) 

a.e. in Q/or alio < t <T. 

Proof: Let io be the supremum of all values of t > for which the statement 
holds, i.e., y(x, t) > z(x, t) a.e. in Q for all < t < to (it is possible that to - 0). 
If to — oo, there is nothing to prove. If to < oo, then by definition 

yMo)>z(-,fo). (6) 

It follows, by Lemma [27T1 that 

C := essinf Q(y(-, to) - z(-, to)) > 0. 

We now define the following sets 

Q+ := {x € Q : y(x, t ) - z(x, t ) > 
Q_ := jx € Q : y(x, to) - z(x, to) < § j . 

By the continuity of the mappings t h» y(-, f) and t h-> z(-, t) from [0, T) to 
L°°, there exists a time interval b\ > such that 

y(x, > z{x, t) for all ^ e [t , t + 6i) for a.e. x e Q+. (7) 

We then turn our attention to the set Q_, where 

ess sup 3>(y(-, t ) - z(-, f )) < ess sup [y(-, f ) - z(-, f )] - 

ess inf Q(y(-, f ) - z(-, f )) < § - C. 
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By the differentiability of the mappings t h-> log y(-, t) and t h-> log z(; t) from 
[0, T) to L°°, there exists for every e>0a time interval 62 > 0, such that for 
all f € [to, t + 5 2 ) , 



||log y(-, f) - log y{-, to) - (t - t )9(a - y(-, £ ))|L < e(t - to) 
||logz(-,f) -logz(-,f ) - (f - t )?(a -z(vfo))|L < e V ~ f o) 
Thus, for t € [to, to + 62)/ 

log ^4 > log - ( f - 'omi/(v to) - z(-, k)) - 2e{t - t ) 

Z\-, t) Z(-, to) 

> -(t - to)ny{; to) - z(; h)) - 2e(t - t ), 

where the last inequality results from ([6]). Choosing e - C/8 and using ^ 
we have that for t 6 [to, to + 62), 

y(-/ 1) c 

Taking 6 = min(<5i, 62) and combining (0 and (O we obtain that 

y(-,t)>z(-,t) for allf e [t 0/ 1 + 5). 

Thus, ([5} holds for all t € [to, to + 6) in contradiction with the definition of 
to, which concludes the proof. ■ 

The comparison principle guarantees the boundedness of y(-, t): 

Proposition 2.3 (Boundedness in L°°) Let y e ^([0, T);L°°' + ) be a solution of 
(TJJ) with initial data yo- Then there exists a constant K > 0, given by (TTT) and 
depending on the initial data, such that 

sup y(; t) <a + Kn(x). (10) 

0<f<T 

Proof: Since ess info n(x) > 0, then there exists, given yo, a constant K > 
such that 

z(x) - a(x) + Kn(x) > yo(x) a.e. in Q. 
Specifically, we can choose 

y (x) - a(x) 

K = ess sup — . (11) 

xen n(x) 

The function z is an equilibrium solution of (fT]), and by the previous propo- 
sition y(-, t)<z for all < t < T. ■ 
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Theorem 2.2 (Global existence) Let yo € L°° ,+ be given. Then (JT) has a unique 
solution y e C\[0, oo), L°°' + ). 

Proof: This is a direct consequence of the short-time existence and unique- 
ness (Theorem 12 and the bound (fTTJb for initial data yo e By the 
continuation theorem for autonomous ODEs, if T < oo and [0, T) is the 
maximal time of existence of the solution y, then 

limsup ||y(-,f)ll°o = oo. 

Since the norm ||y(-, t)\\oo is continuous in time, this violates the bound (flOl) , 
hence the maximal existence time is infinite. ■ 

3 Asymptotic convergence of y(-, t) to M 

Having established the global existence and boundedness of solutions to 
(dJ, we proceed to study the long-term behavior of these dynamics. As 
in the previous section, it is always assumed that system Q} satisfies As- 
sumption [U The first proposition establishes the existence of an integral of 
motion: 

Proposition 3.1 The functional T : L°° /+ — > R defined by 

T(z) := I n(x)logz(x)dy.(x), 
Jq 

is an integral of motion, that is, ify e C 1 (R + ; L°° /+ ) is a solution of CD, then 

r(y(vt)) = r(y ) 

for all t > 0. 

Proof: Differentiating we get 

%-T(y(;t))= f n(x)^^-d^x) = (n,9(a-y(- / t))) = 0, 
at Jq y\x, i) 

where the last equality follows from the symmetry of 7 and the fact that 
n e NCJ>). M 
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The next two propositions reveal the "dissipative" nature of (D through the 
construction of two Lyapunov functionals. Note that by considering the 
equilibrium, y(x) = a{x) + yn(x) for sufficiently large y, we have 

ess inf y(-) > 0, and T(y) = a. 

Proposition 3.2 Let y e C 1 (IR + ; L°°' + ) be a solution of Q} with y(x) defined as 
above. Then the "entropy" functional 

dp.(x) 

is positive and non-increasing in time. 

Proof: The positivity of V a follows from the fact that z - log(z) > 1 for z > 0, 
and the positivity of y(x, t) and y(x). Differentiating along trajectories we 
get 

j-V a [y(; t)] = £ ^| [y(x, t) - y(x)] dp(x) = (y(-, t) - y, 9(a - y(v t))) 
= (?(y(-, f) - y), 7>(a - y{; t))) = -\\7(a ~ y{; < 0, 

where we have used the fact that 7 is an orthogonal projection and 7{y) - a. 
■ 

Proposition 3.3 Let y e C 1 (R + ;L 0O - + ) be a solution of ©. Then the "energy" 
functional 

V h [y(;t)]:=\my(;t)-a)\\ 2 2 . 

is non-increasing in time. 

Proof: By explicit differentiation along trajectories we get 

%V b \y(; t)] = 2 {ny{; t) - a), T(y( v t) 7{a - y{; t)))) 

= -2C?{y{;t)-a),y{;t)?{y{;t)-a)) 

= -2\\y 1/2 (;t)ny(;t)-a)\\ 2 2 <0, 

where we have used the properties of 7 and the positivity of y. ■ 

The identification of the two Lyapunov functionals yields immediately the 
asymptotic convergence of y{-, t) to the equilibrium manifold M. 



V a [y(;t)]:= [ y(x) 
Jn 



log 



y(x, t) 



yfr) 
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Theorem 3.1 Let y e C 1 (R + ;L 0O ' + ) be a solution of ©. Then 

\im?(y(;t)) = a inL 2 . 

t — >00 

Proof: We need to prove that 

(dist L2 (y(-,f),M)) 2 = \my(;t)-a)\\ 2 2 = V b {y{-,t)) 

tends to zero as t — > oo. Since the functionals V a , Vf, are both non-negative, 
bounded from above (Proposiiton I2.3D and non-increasing in time, both 
must converge to limits as t — > oo. Since, furthermore, 

j- t v a [y(;t)] = -v h [y(-M 

the limit of V], must be zero. ■ 



Example: Assume /j(O) = 1 and let 7 be the orthogonal projection in L 2 to 
the space of constants, i.e., 



cp/xx) =/(*)- r f{x')dii( X '), 

Jn 



and a e L°° satisfies 



I fl(x)d l u(x) = 0. 
Jo 

The system Q} takes the form 

j t y(x, t) = y(x, t)(a(x) + y(x', t) d^x') - y(x, t)), (12) 

with initial condition y(-,0) = J/o e L°°' + . Theorem l2.2l asserts the existence of 
a global solution y e C 1 (R + ; L°°' + ). By Proposition |2.3l there exists a constant 
K > such that 

sup y(-, < a + K. 

f>0 

Finally, by Theorem |3.1[ 

lim(y(-, t) - I y{x' ,t)dn(x')\ - a inL 2 . 
^°°\ Jo / 
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4 Asymptotic convergence of y(-, t) 



We now question under what conditions does y(-, f) converge, as t — > oo, to 
a specific equilibrium in M. Note that the L 2 -convergence of y(-, f) can be 
decomposed into 

lim y(-, f) = lim :P(y(-, t)) + lim Q(y(-, f)), 

t—*oo t — >oo t—*oo 

where 

Q(y(-,0) = (n,y(v,f))n. 

We have just proved that the first term on the right-hand side converges to 
a. It remains to verify under what conditions 

£(*):= ((y(-,t),n) (13) 

converges as f — > oo. 

Since, on the one hand, M consists of functions of the form a{x) + a n(x), 
for some a e R, and on the other hand, by Proposition 13.11 the functional 
T(y(-, t)) is conserved, the existence of a limiting solution in M requires the 
following assumption: 

Assumption 2 There exists some y* 6 M such that 

I n(x) logy Q (x)dn(x) = I n(x)logy*(x)dp(x). (14) 
Jo Jn 

Assumption |2] is a restriction on the initial conditions yo- It assumes the 
existence of a constant a which solves the equation 

I n(x) log [fl(x) + an(x)] dy.{x) = I n(x) log y (x) d/.((x), (15) 
Jn Jn 

under the constraint that ess info[fl(x) + an(x)] > 0. 

If we define the set C c IR by 

6 = U 6 R : ess inf [a(x) + £n(x)] > o) (16) 
I xen J 

and O : C — > R by 

®(Q= [ n{x)log[a{x) + Zn{x)]dp{x), (17) 
Jn 
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then Assumption[2]is equivalent to the statement 

n(x) logyo(ac) dp.{x) e 0(6). 

Note that S is in fact an unbounded interval, for £ 6 C implies that ^ e C 
for all & > £. 

The next proposition shows that such an a, if it exists, is unique. 

Proposition 4.1 Given an initial data yo e L°°' + , the function y* satisfying As- 
sumption^ if it exists, is unique. 

Proof: Uniqueness follows at once from the fact that 
for all a e Q. ■ 

Example: Consider again the example from the previous section. For con- 
creteness set Q = [0,1], with y. the Lebesgue measure and a(x) = sin27ix. 
Then, since n = 1, the equilibria in M consist of functions of the form 

sin 2nx + a, 

where a > 1, i.e., C = (1, oo). For a € C, 

O(a) = I log[sin 2nx + a] dx > - log 2. 
Jo 

It follows that Assumption|2]is satisfied if and only if 

l 

\ogy (x)dx > - log 2. 



The following proposition asserts that the convergence of y(-, t) is guaran- 
teed if the solution remains bounded away from the boundaries of the cone 
of positive solutions L°° ,+ . 



I 



I 
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Proposition 4.2 If 



liminf essinfy(x, t) > (18) 

f->oo xeCL 



then Assumption's satisfied. Moreover, 

lim |3(£) = a, 



t — >oo 



where p(t) is given by ((131 ) and a is the (unique) solution to ((151 ). Thus, y(-, t) — » y* 
m L 2 , zy/zere y* = a + an. 

Proof: Take any sequence of times t m that is increasing to infinity. Since 
y{-,t) is uniformly bounded in L°° (Proposition 12.31 ), then |3(£) is bounded, 
and there exists a subsequence t m such that fi{t m ) converges to a limit y, 
hence 

Urn Q(y(-,f m J) =yn inL°°. 
Theorem 13.11 implies that 

lim \ y(; t mk ) - Q(y( v t m J) - a] = in L 2 , 
from which follows that 

lim y(x,t mk ) = a(x) + yn(x) 

J;— >oo 

in L 2 , and so it has a sub-subsequence y{-,t mk ) which converges a.e. in Q. 
Note that ( flBl implies that a.e. a(x) + y n(x) > 0. This implies that 

lim n(x) log y(x, t m ) = n(x) log [a(x) + yn(x)] 

j—>co J 

a.e. Moreover, from ( fT8|) and the fact that y(-, £) is uniformly bounded we 
also have 

sup||logy(-,f)IU < oo. 

f>0 

Using Lebesgue's dominated convergence theorem we conclude that 

lim I n(x)\ogy(x,t mk ,)d[i(x) = I n(x)log[fl(x) + yn(x)] dy.{x). 
H°° Ja. ' Ja. 

By Proposition ^. 1 1 we have for all m, 

I n(x)logy(x,t m )dy.(x) = I n(x)logy (x)d[i(x), 
Jn Jn 
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therefore 

I n(x) log [a(x) + yn(x)] dp{x) = I n{x)logy (x)dy,{x). 
Jq Jn 

Thus, Assumption [2] is satisfied and it follows, by the uniqueness of y* , 
hence the uniqueness of a in ((15)) , that y = a. We have shown that every 
sequence f}(t m ) has a subsequence j3(i W)fc .) which converges to a. It follows 
from an elementary theorem of calculus that fi(t) tends to a as t — > oo. This 
completes the proof. ■ 

Note the immediate corollary: 

Corollary 4.1 If Assumption\2\does not hold then 

lim inf ess inf y(x, t) - 0. 

f-»oo xeO 



Condition ((18)) is a sufficient condition for y(-, to asymptotically converge 
to an element of M. The problem is that it is a property of the solution, and 
it is not clear a priori when does it hold. In the remaining part of this section 
we establish two situations for which ( fT8l holds. In the first case yo has to 
be sufficiently large in the following sense: 

Proposition 4.3 If there exists a constant K such that 

yo(x) > a(x) + Kn(x) > a.e. in Q, 

then condition dT8t holds. 



Proof: This is an immediate consequence of the fact that a + Kn is a stationary 
solution of (d), and the comparison principle (Proposition l2.2|) . ■ 

The second situation that can be analyzed is when a and n are simple 
functions, i.e., they have the form 

m m 

«(*) = Yu a iXnt( x l n ( x ) = "<<W x )' 
i=l i=l 

where Qi, . . . , Q m is a measurable disjoint partition of Q. 
Proposition 4.4 If a and n are simple functions then ((18)) holds. 
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Comment: The implication of this proposition is that (18} holds for any 
finite-dimensional approximation of <(XJ) . In particular, the solutions to dis- 
crete approximations of (D with positive initial data always tend to equi- 
librium solutions as t — > oo. 



Proof: We first prove the proposition for the particular case in which yo(x) = 
c > (a constant function). Note that if yo and n are simple functions with 
respect to the partition (Q;), then the right hand side of ((T]) is also a simple 
function, in which case y(x, t) is a simple function, constant on each of the 
sets Q,, for all t > 0. We denote by y,(£) the restriction of y(x, t) to the set Q,. 

Let M be a bound on \y{x, t)\ (such a bound is guaranteed to exist by Propo- 
sition Then for all t > 0, 



f 

Jn 



n(x) log y(x, t) dp(x) = ^(Q;) m log y f (t) 

!=1 

< log M I n(x) d^t(x) + I min my.{D.i) I log I inf y(x, i) 

Jq \l<i<m I \xeCl ) 

On the other hand, by Proposition 13. II 

I n(x) log y(x,t)dy.(x) = I n(x)logy (x)^(z) = logc I n{x)d[i{x), 
Jo. Jn Jq 



hence 

(log c - logM) J Q n(x) djj(x) 



inf y(-, t) > exp 



mini< ! < m [n,fi(Q/)] 



> 0. 



This completes the proof in the case of constant initial conditions. The 
general case follows at once from the comparison principle, as any solution 
with initial data yo € L°° ,+ can be bounded from below by the solution for 
constant initial data c — ess infQ yoM- * 



5 Discussion 

We studied a class of quadratic evolution equations, inspired by models of 
viscoelastic fluids. Motivated by the physical model, we considered initial 
data in the cone of positive functions. We showed that the cone of positive 
L°° functions is an invariant set, and that solutions in this set exist for all 
times. As t — > oo the solutions tend, in the L 2 -norm, to the equilibrium 
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manifold M. The convergence of solutions to specific equilibria in M could, 
however, only be proved under additional assumptions. 

The following points remain open: (i) Do solutions always tend to a specific 
equilibrium if Assumption^ is satisfied? We were unable to prove it, nor to 
find a counter example, (ii) Do solutions converge, as t — > oo, in situations 
where Assumption [2] does not hold? While, in such case, the solution 
cannot converge to an equilibrium in M (Corollary 14.1) 1. it can, in principle, 
converge to an equilibrium on the boundary of the cone, 

L^ = {yeL°°: y(x)>0}. 

(iii) Does the solution converge to M in any LP -norm, for p > 2, and in 
particular, for p — oo? 

Another question is whether our results remain valid when the kernel of 
the projection T has dimension greater than one. The comparison principle 
(Proposition [221) no longer holds in this case, and as a result, we no longer 
have a bound on the L°° norm, nor do we have a global existence theorem. 
Assuming, however, that a solution does exist for all times, it is easy to see 
that Proposition |3.3l still holds, i.e., the "energy" functional is a Lyapunov 
functional. To prove that the "entropy" functional V a is also a Lyapunov 
functional, we need to have a positive function y such that !P(y) = a. If such 
function exists then Proposition |3.2| remains valid, and T(y) tends to a in the 
L 2 -norm (Theorem 13 .!)) ■ 

System (TJJ) can be generalized in many different ways, for example, with 
y being a matrix valued function and products reinterpreted as matrix 
products; this is indeed the appropriate setting in the viscoelastic context 
PI. Another generalization of Q]) is when T is a general non-negative 
operator (not necessarily a projection), i.e., (y, T(y)) > for all y € L 2 . We 
believe that such a system still exhibits global-in-time existence for positive 
initial data, as well as asymptotic convergence. 
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